Scattering of quantum wave packets by shallow potential islands: 

A quantum lens 
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We consider the problem of quantum scattering of a localized wave packet by a weak Gaussian 
potential in two spatial dimensions. We show that, under certain conditions, this problem bears 
close analogy with that of focusing (or defocusing) of light rays by a thin optical lens: Quantum 
interference between straight paths yields the same lens equation as for refracted rays in classical 
optics. 

PACS numbers: 03.65.Nk, 03.65.Sq, 42.25.Fx 



I. INTRODUCTION 

The intrinsic connection between the motion of classi- 
cal particles on one hand and the propagation of quantum 
matter waves on the other has occupied minds of scien- 
tists since the early days of quantum theory. De Broglie 
was one of the first to realize that "for both matter and 
radiations ... it is necessary to introduce the corpuscle 
concept and the wave concept at the same time" 
Subsequently, invaluable contributions of Ehrenfest, Van 
Vleck, Feynman, Gutzwiller, Maslov, among many oth- 
ers, shaped our current understanding of quantum wave 
propagation in terms of interference of classical trajec- 
tories. However, a number of important questions con- 
cerning quantum-classical correspondence remain open. 
These questions fall under the scope of the area of math- 
ematical physics known as quantum chaos 043- 

The motion of a classical particle can be conveniently 
described by means of a phase-space trajectory. The 
Heisenberg's uncertainty principle however does not al- 
low for the notion of the classical trajectory to be directly 
carried over to quantum theory: the particle's position 
and momentum can not be specified simultaneously. 

One natural extension of the classical concept of a 
point in the phase space is provided by a localized quan- 
tum wave packet that can be parametrized by its mean 
position and momentum, and dispersion quantifying the 
phase-space extent of the wave packet. According to the 
Ehrenfest theorem the time evolution of the mean 
position and momentum is governed, for short enough 
times, by the classical equations of motion. In other 
words, the wave packet center follows the corresponding 
classical trajectory. 

An issue of the wave packet spreading, i.e., how the dis- 
persion depends on time, is however much more complex. 
Loosely speaking, there are two main mechanisms of the 
spreading: (i) a classical-like broadening of the wave 
packet due to forces exerted by an external potential, 
and (ii) an intrinsically quantum-mechanical spreading 
dictated by the uncertainty principle. Due to the inter- 



ference nature of quantum dynamics, the overall spread- 
ing is not a simple "sum" of the two contributions, but 
rather a more intricate process. 

A natural question arises: is there an intuitive and, at 
the same time, quantitative theoretical description of the 
phenomenon of quantum spreading? In this paper, we de- 
velop such a description, based on the short-wavelength 
approximation to quantum dynamics, for the simple sys- 
tem of a two-dimensional quantum wave packet scattered 
by a weak Gaussian potential. In particular, we show 
that the quantum scattering process bears a close math- 
ematical analogy with the phenomenon of focusing (or 
defocusing) of light rays by a thin lens, and can be de- 
scribed using the "language" of geometrical optics. In- 
terestingly, on the quantum side the use of the Eikonal 
approximation Q, i.e. including interference of straight 
paths, yields the same thin lens equation as derived in 
classical optics from refracted light rays. 

Our theoretical approach provides an intuitive picture 
of the wave packet spreading 0, Q, and its quantita- 
tive predictions are found in good agreement with results 
of an "exact" numerical solution of the time-dependent 
Schrodinger equation. 



II. THEORY 

We consider a quantum particle of mass m that evolves 
in the two-dimensional position space under the influence 
of the external potential 



V(q) = V Q e 



q Aq 



(1) 



Here, Vq quantifies the strength of the potential, q is a 
column vector representing the particle's position, and 
A is a 2-by-2 orthogonal matrix with eigenvalues a\ and 
02 corresponding, respectively, to orthonormal (column) 
eigenvectors ei and e 2 . In other words, 



A = (ei e 2 ) diag(ai,a 2 ) (ei e 2 ) 



(2) 
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where |ei| = |e2 1 = 1 and ei -e2 = 0. Hereinafter, the dot 
"•" stands for the scalar product, and the superscript " T " 
denotes the matrix transposition. Equation (p} describes 
a Gaussian potential "island" centered at q = 0, see 
Fig. Q] The spatial extent of the island is characterized 
by the length 1% = in the direction of the vector 

ei and by li — l/y/a~2 in the direction of &2- 




ei 



FIG. 1: (Color online) Schematic illustration of the scattering 
system under consideration. 

The dynamics of the quantum particle is fully de- 
scribed by its time-dependent wave function, ^t(q). The 
latter is related to the initial wave function, ^o(q), by 
means of a propagator K in accordance with 



= J dq'^(q,q';t)vl> (q'). 



(3) 



Here, the q'-intcgration runs over the whole two- 
dimensional plane. In general, the propagator K is deter- 
mined by solving the time-dependent Schrodinger equa- 
tion, 



F(q) tf(q,q';*) = 0, (4) 



subject to appropriate boundary and initial conditions, 
e.g., see Ref. 9]. Generally, this task is analytically 
formidable, and one is left to resort to numerical com- 
putations. In certain cases however analytical progress is 
made possible by constructing "sensible" approximations 
to the propagator. In this paper, we focus on one such 
approximation commonly referred to as semiclassics. 

The semiclassical (or short-wavelength) approximation 
of the propagator K(q, q';t) is formulated in terms of 
classical trajectories 7 that start at the point q' at 
time and end at q at time t. More precisely, 7 = 
•jr(r), r € [0,t]}, such that + VV(r) = with 

r(0) = q' and r(t) = q. The approximate propagator, 
commonly referred to as the Van Vleck-Gutzwiller prop- 
agator, can be written as 



2nih 



exp 



1 in 
-5 7 - Vl 



(5) 



where 



S 7 = S[r](t) 



/^{f(fw) 2 -^» 



(6) 



is Hamilton's principle function along the trajectory 7, 

D 1 = |det(-V q <V q S 7 )| (7) 

is the stability factor of 7, and v 1 is the so-called Maslov 
index, counting the number of conjugate points along 7; 
as far as our problem is concerned, the Maslov index is 
identically zero, = 0. 

The approximation ([5]) is known to be reliable when 
applied to wave functions ^(q) representing a quantum 
particle of sufficiently high kinetic energy Eq. We ad- 
ditionally assume that the Gaussian potential, given by 
Eq. ([T]), is weak compared to the particle's kinetic en- 
ergy, I Vo| <C E . The validity of this condition is at 
the heart of our "shallow potential island" approximation 
corresponding to the Eikonal approximation in scatter- 
ing theory [6]. In this case (see Appendix |A"|) . 5 7 can be 
well approximated by S l01 where 70 is the straight, free- 
particle trajectory leading from q' to q in time t, i.e., 
7o = {qr/t + q'(l - r/t) , r e [0, t]}. Thus, we write 



c m i 



,'|2 



-| t drF(q' + ^(q-q') 



(8) 

The integral in the right-hand side of Eq. ([5} can be 
straightforwardly evaluated to equal 



2A 



■ exp 



(q- Aq)(q'-Aq')-(q-Aq') : 



A 2 



erf 



q • A(q-q / ) 

A 



— erf 



q , A(q -qQ 
A 



with 



A = ^(q - q') • A(q - q') . 



(9) 



(10) 



Expression ©, and therefore Eq. ©, can be further sim- 
plified by taking into account the identity 

(q-Aq)(q'.Aq')-(q-Aq') 2 = |qxq'| 2 detA, (11) 
where "x" denotes the vector product. This yields 



|q- q 



erf 



2V 1 ~* 1 2A 
q • A(q - q') 



A 



exp 
erf 



|q x q'| 2 det A 
J? 

q'-A(q-q') 



A 



(12) 



A substitution of Eq. (fT2"|) into Eq. ([SJ leads to an ex- 
plicit, closed-form expression for the semiclassical prop- 
agator, _K"(q, q',t), and, therefore, provides the complete 
solution of the time-dependent scattering problem in the 
short- wavelength regime. 
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The expression for the propagator, K(q, q', f), becomes 
especially simple and allows for an intuitive interpreta- 
tion in the following special case. Let us consider a setup, 
in which the "receiver" q and the "source" q' lie on the 
opposite sides of and at almost the same distance, large 
compared to h, from the center of the Gaussian scatter- 
ing potential, and in which the vector (q — q') is nearly 
aligned with one of the principal directions (taken, for 
concreteness, to be ei) of the potential island. In other 
words, we are interested in the asymptotic form of the 
function K(q, q',i) in the case that 



Lex + £ • 



and 



|£|,|£'U«£- 



(13) 



(14) 



Substituting Eq. ([13) into Eqs. (JTOJ) and ([12]), taking 
into account that erf (±z) —¥ ±1 as z — > +oo, and keeping 
only terms to the leading order in \£\/L and |£ \/L in the 
argument of the exponential function, we obtain 



S 7 ^|2Le 1+ £-^| 2 

-—L VoteXP { W 



(15) 



Then, using the basis representations £ = Ci e i+?2 e 2 and 
£' = ^ei + £ 2 e 2 in Eq. (fT5|) . and further assuming that 



we write, approximately, 



S " , 



2 L 



v t + + 



where 



(16) 



(17) 



(18) 



and 



= - Q 2 + ^r^Vot (6 + Q 2 ■ . i - 



8 L(l 2) 



In view of Eqs. (fl~7HT9l) . the stability factor along 

^)iy Dry Dry 



the trajectory 7 can be written as D 1 ~ DjDi?\ with 



D« = = m/t and Z^ 2) = |-%%sf } | = 

m/t - v^iV t/(4L(/ 2 ) 2 ), provided that 



t 2 < 



(20) 



I Vo| Zi 

In fact, Eq. ([20]) states a necessary condition for Eq. ([T9j) 
to constitute a "healthy" perturbative expansion of S\ 
in powers of Vq. 



The full semiclassical propagator, Eq. ([5]), can now be 
written as 

X(q,q,i)^exp T — 

x tf (£ + &,-£ + &*) #v(6,&t), (21) 

where 

#0(2, 2', t) 



J^expA — (z-z') 2 ) (22) 

is the free-particle propagator describing the motion of 
the particle in the ei-direction, while 



K v (z,z',t) =K (z,z',t)Ji 



yfth Vqt 2 



x exp 



v 4 Lm{l 2 ) 2 

nh V t {z + z' f 



8 L h 



(h) 2 



(23) 



accounts for the wave function spreading in the orthogo- 
nal, e2-direction. Clearly, Ky — > Ko as Vb — > 0, recover- 
ing the free-particle limit. 
We now observe that 



+00 



K v (z,z',t) = / d£K (z,{,T/2) 

J —OO 

^hV T / J^h V T 2 



0X1> " 2 L h 



1 



4 L m(/a) 2 



(24) 



Equation (1241) is an identity, and can be verified straight- 
forwardly by evaluating the Gaussian integral in the 
right-hand side. Then, after Eq. (|20[) is taken into ac- 
count, Eq. ([2~4"]) reduces to 



K v (z,z',t) 



d(K (z,(;,T/2) 



x exp 



.v^hVor e 



2 L H (25) 

The physical picture offered by Eq. (I2"5l is as follows. 
The propagator Ky, evolving a quantum state during 
time t, can be view as a result of three consecutive opera- 
tions: (i) a free-particle propagation during time r/2, (ii) 
an instantaneous phase change, or "kick" , of the quantum 
state, and (hi) another free-particle propagation during 
r/2. This interpretation, and the physical meaning of 
the kick operator, becomes apparent when the evolving 
quantum state is given by a Gaussian wave packet. To 
this end, we consider as initial state the two-dimensional 
wave packet 



M&&) = 4 1) (&Pl)4 2 \&P2) (26) 



with 



1 

1 

ira 2 



exp 



.mv 
mv z 2 \ 



+ z 



(27) 
(28) 
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Here, pi and p 2 (usually termed radii of curvature 
are two, generally complex-valued, parameters. The real- 
valued function a — cr{p), defined in accordance with 



1 



mv ?sp 



ir, = ^ ~ =~ i^, (29) 



quantifies the position-space dispersion of the wave 
packet. Furthermore, v specifies the average velocity 
(and mv the average momentum) of the particle. 

Now, acting with the propagator (l2"TT) on the initial 
state (j2l)|) , which is assumed to be spatially localized 
around the position vector q' = —Lei (or around the 
origin in the ^'-coordinate frame), we obtain the quan- 
tum state after time t locally, in the vicinity of the point 
q = Lei (or around the origin in the ^-coordinate frame): 



where 

^\z;p) 



■ y^F H Vnt 



*^4 1 '(^Pi)^ 1 (^,P2), (30) 



(2), 



dCKo(L + z,-L + C,t)tl>g\<;;p), 



^ (2) (z;p) S / dCK v (zX,t)4 2 \C;p). 



(31) 
(32) 



As we are concerned with the semiclassical limit, it is rea- 
sonable to expect that, in the course of its time evolution, 
the quantum wave packet remains concentrated around 
the corresponding classical trajectory. This means that 
the center of the wave packet, starting from the point 
—Lei at time 0, reaches the point Lei hi time t, such 
that 



vt = 2L. 



(33) 



It is at this instant that ^4(^1,^2) is localized around the 
origin in the ^-coordinate frame, and that the propagator 
approximation, given by Eq. (|2f \ , proves the most useful. 

Fixing the time t in accordance with Eq. (|33p and eval- 
uating the Gaussian integral in Eq. (I3f [) . we obtain 



(34) 



where fa — E t/h — (1/2) arg(l + vt/pi), and 



En 



(35) 



denoting the kinetic energy of the corresponding clas- 
sical particle. The physical interpretation of Eq. (|3"4")l 
is that the ei-component of the wave packet retains its 
Gaussian shape, as its center travels in space on top of 
the corresponding classical trajectory (and in agreement 
with the Ehrenfest theorem). The spreading of the ei- 
componcnt of the wave packet is entirely described by 
the linear transformation of the corresponding radius of 
curvature, pi — > pi +vt, and, in the weak potential limit, 
this spreading is not affected by the external potential. 



We now focus on the time evolution of the e 2 - 
componcnt of the wave packet. As before, we keep the 
time t fixed in accordance with Eq. (|33p . Substituting 
Eq. (1231) into Eq. (|32l) . and successively evaluating two 
Gaussian integrals, we obtain 



4 2 \z;p 2 ) = e^ ^ 2) (z;p 2 ), 



where p' 2 



■vt/2, 
1 

P+ 



p- = P2 + vt/2, 



f 



1 1 

- + 1 

1 E {l 2 f 

V^Vq h 



(36) 



(37) 



(38) 



and cj) 2 = -(1/2) [ arg(p 2 /p + ) + arg(p_/p 2 )] . 

The physical picture of the wave packet spreading, 
offered by the central Eqs. (15fJH3"8"|) . bears close anal- 
ogy with the focusing (or defocusing) of light rays 
by a thin optical lens. Indeed, the well-known thin 
lens equation, [object distance] -1 + [image distance] -1 — 
[focal length] -1 , can be readily recovered from Eq. ([37jl 
by interpreting — p_ and p+ as the "distances" from, re- 
spectively, the object and its image to a lens of the focal 
length /; the role of the lens is played here by the Gaus- 
sian potential island, see Fig. [TJ 

It is important to point out that the above analogy 
between the wave packet scattering in quantum mechan- 
ics and the ray focusing in optics is not a trivial one: in 
the quantum-mechanical case, the "distances" — p_ and 
p + are intrinsically complex-valued and can only be re- 
lated to the true distances, encountered in optics, in a 
nonlinear way. 

Finally, we note that our simple wave-packet- 
propagation construction can be generalized to arbitrary 
times t. This generalization can be summarized as fol- 
lows. Suppose that the initial wave packet is centered 
around a point with coordinates (Q,0), where Q < 0, 
and that the initial wave function is given by the prod- 
uct ip(p(qi — Q; pi) V'o 2 ' (92; P2)- Then, at a later time t, 
the wave function is given, up to an overall phase factor, 
by ipo\(li - Q-vt;pi + vt)4 2) {q 2 ; p' 2 ), where 



P2 



P2 

f>+ 



vt , 
-v(t 



\Q\M 



t<\Q\/v 
t>\Q\/v 



(39) 



and p+ is determined from Eq. (|37|) with p_ = p 2 



III. NUMERICAL CONFIRMATION 

We now confirm the validity of the "quantum lens" 
formulae, given by Eqs. (|57|) . (l3"5l) . and (|31?j) . by com- 
paring their predictions to results of numerical simula- 
tions. The latter were performed by solving the time- 
dependent Schrodinger equation, governing the evolu- 
tion of the wave packet, numerically, using the method 



5 



of expanding the propagator in a series of Chebyshev 
polynomials of the Hamiltonian. The reader is referred 
to Refs. [IiTll3l ] for a comprehensive description of the 
method and its implementations. 

The concrete system that we consider is schemati- 
cally illustrated in Fig. Q] and described by the following 
set of parameters. Hereinafter, we adopt atomic units, 
ft, = m = 1. The external Gaussian potential, Eq. (fl}, 
is characterized by = 0.1 and l 2 — 1, and the po- 
tential strength Vq plays the role of a variable param- 
eter, with values ranging between 10 and 40. The ini- 
tial state of the particle is given by the wave function 

wHqi - Q',Pi) ^ < £\<12;P2) with Q = -0.8, v = 60, and 
Pi = p 2 = —i(mv J 'H)(ao) 2 , where er = 0.1 quantifies the 
initial position-space dispersion of the wave packet, cf. 
Eq. (|2U)) . Note that the kinetic energy of the classical 
particle Eq = mv 2 /2 = 1800 is large compared to the 
strength of the external potential. 



,x 1CT 



calculation yields 
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FIG. 2: (Color online) Spreading of the wave packet in the 
direction orthogonal to the direction of propagation. See text 
for details. 



Our aim is to study the wave packet spreading in the 
direction orthogonal to the propagation direction. This 
spreading is given by the time-dependence of the disper- 
sion along the e2-axis, 




dq(<? 2 ) 2 |* t (q) 



(40) 



(Note that, due to the symmetry of \&t(q) under the re- 
flection q 2 — > — <72, the expectation value of <j2 is zero, 
i.e., / dq<72|^t(q)| 2 = 0.) For our choice of the initial 
state, (72 = &o at t = 0, and oi increases as a function 
of time. In the limit Vq = 0, this increase is determined 
by the free-particle spreading, cr| rec = <r{p 2 + vt), with 
the function a(p) defined by Eq. (I29p . A straightforward 



free 



\mao _ 



(41) 



In the case Vq ^ 0, the external potential causes addi- 
tional spreading that can be quantified by 



Acr 2 = ct 2 



free 



(42) 



Figure [2] shows the dependence of A&2 on time t for 
four different potential strengths, Vq — 10 (blue), 20 
(magenta), 30 (green), and 40 (red). The solid curves 
represent the results of the "exact" numerical solution 
of the time-dependent Schrodinger equation, while the 
dashed curves show the analytical "lens" approximation, 
namely <j(p' 2 ) — oi' 00 with p' 2 , calculated in accordance 



with Eqs. ([M]), ((33, and J3S]). As expected, Aa 2 ^ 
for t < \Q\/v ~ 0.013, corresponding to the time that 
it takes for the classical particle to reach the potential 
island. Figure [5] shows the theoretical predictions to be 
in a reasonable agreement with the numerical results. It 
also confirms that the agreement improves as the poten- 
tial strength is decreased. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper, we have constructed an approximate an- 
alytical solution to the problem of a scattering of a lo- 
calized quantum wave packet by a weak Gaussian poten- 
tial. Our solution is valid in the semiclassical regime, in 
which the particle's de Broglie wavelength can be con- 
sidered short compared to all other length scales of the 
system. We have shown that the quantum scattering 
process is closely analogous to the phenomenon of focus- 
ing (or defocusing) of light rays by a thin optical lens. 
In particular, the mathematical formula quantifying the 
wave packet spreading, Eq. (|37[) . is largely equivalent to 
the thin lens formula of geometrical optics. The main 
difference between the thin lens formula in optics and 
Eq. (f3"T)) is that the former operates with true real-valued 
distances from the lens to the object and to the image, 
while the "distances" p_ and p + entering Eq. (|3T|) are 
intrinsically complex- valued. It is only in the classical 
limit, mv/H — ¥ oo, that p_ and p + become real- valued, 
and that the optical thin lens formula is recovered. 

It is instructive to further compare the ray optics pic- 
ture with our quantum mechanical result. The thin lens 
formula in optics, which is the classical limit of Eq. (|37l) . 
effectively describes the deflection, or bending, of light 
rays (classical trajectories) , induced by the lens (external 
potential). However, only straight, unbent trajectories 
have been used in our semiclassical derivation of Eq. ([3T]) . 
This seeming paradox is resolved by the following ar- 
gument, originally presented in Ref. 14| and for read- 
ers' convenience reproduced in Appendix [A] The main 
building block of the semiclassical propagator, Eq. ([5]), 



6 



is the Hamilton's principal function along the classical 
trajectory connecting the initial and final points of the 
propagation. However, in a sufficiently weak external 
potential, the value of the Hamilton's principal function 
along the true (generally bent) classical trajectory is very 
close to that along the corresponding straight (unbent) 
trajectory, making the precise geometrical shape of the 
trajectory unsubstantial. 

The approach taken in this paper is conceptually sim- 
ilar to the one used in Ref. [1(| to analyze the spreading 
of quantum wave packets in the Lorentz gas. The latter 
consists of a particle moving in an array of fixed elastic 
scatterers, taken to be hard disks (spheres) in two (three) 
spatial dimensions. However, it is important to point out 
that in the Lorentz gas, unlike in the system addressed in 
the present paper, one must take into account deflections 
of classical trajectories in order to obtain the quantum- 
mechanical equivalent of the circular (spherical) mirror 
formula. 
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Appendix A: Expansion of Hamilton's principal 
function 

The following discussion is based on the argument that 
was, e.g., presented in Ref. 



Let us consider the Lagrangian 

A[r](r) = |^(r)) 2 - e U(r(r)) (Al) 

along a trajectory r(r). Here e <C 1 serves as a dimen- 
sionless strength of the potential. The corresponding 
Hamilton's principal function is given by 

t 

S e [T](t) = Jdr£ e [T]{r). (A2) 
o 

We now denote by r e (r) a trajectory that satisfies the 
boundary conditions r e (0) = q' and r c (t) = q, and makes 
the action S e [r](f) stationary, i.e., 

^[r e ](r) = 0. (A3) 
or 

Then, for a trajectory i"o(t), that satisfies the same 
boundary conditions, r (0) = q' and r (t) = q, and that 
is the stationary trajectory of 5o[r](t), we have 

S e [r ](t) = S e [T e ](t) + [ dr &r e ] • (r - r e ) 
Jo or 

t 

+ 1 |dr(r -r e ).^[r £ ](ro-r e ) + .... (A4) 
o 

Substituting Eq. (|A3j) into (|A4p and taking into account 
r £ = ro + 0(e) we obtain 

S e [r e ](t) = S e {v Q ](t) + 0(e 2 ) . (A5) 
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